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Abstract 



> 



O 

\^ • We study the transverse momentum distribution of vector mesons produced in uhraperipheral 

o ■ 

■ relativistic heavy ion cohisions (UPCs). In UPCs there is no strong interaction between the nuclei 

o _ 

and the vector mesons are produced in photon-nucleus collisions where the (quasireal) photon is 
emitted from the other nucleus. Exchanging the role of both ions leads to interference effects. 
^ • A detailed study of the transverse momentum distribution which is determined by the transverse 

momentum of the emitted photon, the production process on the target and the interference effect is 
^ ' done. We study the total unrestricted cross section and those, where an additional electromagnetic 

excitation of one or both of the ions takes place in addition to the vector meson production, in the 
latter case small impact parameters are emphasized. 

PACS numbers: 



p. 
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Due to the strong electromagnetic fields surrounding the heavy ions in relativistic colli- 
sions, RHIC and LHC can be seen as a factory of quasireal photons of high energies. One of 
the interesting photonuclear processes studied in these "ultraperipheral collisions" (UPC) 
is the coherent production of vector mesons, in particular p°, which has been measured 
recently at RHIC . The coherent production was identified through the transverse mo- 
mentum distribution of the meson, which is enhanced for values of the transverse momentum 



v± 



< 



1/R where R denotes the nuclear radius. We give a careful theoretical study of the 



process 



(1) 



with ("y4*") and without ("A") an electromagnetic (GDR) excitation of either one or both 
ions. This is of interest for the analysis of the RHIC experiments, as well as, for future 
experiments at LHC where also other vector mesons like {J/"^ and even T) can be studied 
^, y, 0]. While the theory of UPC is generally in a good shape 0, Q, 0], the specific 
question of the transverse momentum distribution has been paid attention to only in less 
rigor 0, l3, 11|. 

Heavy ion scattering offers a unique possibility to study an important interference effect 
[lol | . As will be shown below, the transverse momentum distribution is very sensitive to this 
effect. It is the purpose of this letter to give a careful study of this transverse momentum 
distribution. 





(a) 



(b) 



FIG. 1: A schematic Feynman diagram for the vector meson production in ultraperipheral heavy 
ion collisions (a). The corresponding exchange diagram is also shown (b). 



The kinematics of the process given in Eq. (PJ) is denoted by (see Fig. P) 

p + k^p' + k' + v. 



(2) 



Due to the additional elastic photon exchanges which are schematically denoted by the 
open blobs in Fig. ^ the momenta Q and A are not related to the asymptotic momenta 



by Q = p — p' and A = k — k' as, e.g., in the pp case (without rescattering) For 
small transverse momenta the longitudinal components of the photon momentum and the 
momentum transfer from the vector meson production ("Pomeron momentum") are given 
in the cm. system by the mass my and the rapidity Y of the produced meson as 

Qo = PQ. = Ao = = ^e-\ (3) 

The momenta p and k, see Fig. Q^a), are given by p = m^'u+ (ion 1) and k = rriAU- (ion 2), 
where uia is the ion mass and u± = 7(1, 0, 0, ±/?). In the exchanged process, see Fig. Wljo), 
the photon is emitted from ion 2, the "Pomeron" from ion 1. 

Due to the large value of the Coulomb parameter ri = ^^fi'^^ we can use the semiclassical 

unn 

approximation |2, llj, ll^. We also show how this can be derived from eikonal/Glauber 
theory, more details will be given in a forthcoming publication Using a simple model 
for the meson production process we are able to give analytical results. Implications for 
the current experiments at STAR/RHIC and for future experiments at the LHC jj, 10 1 are 
discussed. In the semiclassical approximation the two ions move along a straight line and the 
process is described by an impact parameter dependent amplitude a{b,v±,Y). In contrast 
to the momentum of the vector meson the momenta of the outgoing ions are not detected 
and the differential cross section is given by 



Irr^Y^f d'bMb,v^,Y)\\ (4) 



d^v±dY 2(2 

The integration over the impact parameter b corresponds to an integration over the unob- 
served momenta k' and p' of the scattered ions. The different processes which can occur 
according to Eq. (P) factorise jl^ : 

afi{b,v±,Y) = anuci{b) ai{b) 02(6) av{b,vuY). (5) 

The strong absorption due to the interaction of the ions for b < 2R is given by anud{b) ~ 
0(6 — 2R) with the nuclear radius R = 7fm. av{b,v±,Y) describes the vector meson pro- 
duction. Additional electromagnetic excitation amplitudes of ion 1 and/or 2 are denoted by 
ai{b) and 02(6). What is chosen for them depends on the additional triggering condition 
used in the experiment. The cross section can be written as 

d^a 1 



dW 2(27r)3 J,^ 



poo 

/ d%Ub)J2Wvib,v,Y)\\ (6) 

J2R 



where fij{b) takes the triggering condition of the measurement into account. In the case 
where the process is observed without any further condition imposed on the excitation of the 
ion(s), we have foo{b) = 1. If the vector meson production together with the electromagnetic 
excitation is measured one uses either /io(&) = -Pi(^) = or /oi(&) = P2{b) = |a.2(&)P 

if the excitation of one of the ions, /ii(&) = Pi{h)P2{h) if the mutual excitation of both ions 
is triggered on. 

The electromagnetic excitation probabilities are given by Pi{h) = with S ~ 5.45 x 
10-^ Z^NA-^/Hm^ [16]. This equation is vahd for 2R < b < bmax = but this cutoff 

can be safely neglected here, since the vector meson production probability falls off at least 
as fast as 1/6^. 

n the semiclassical treatment the amplitude of an electromagnetic process can be written 

= / ^,AUb,Q)MQ), (7) 

AUb,Q) = 2nZeQ'i6{Qu+)-^^^exp{-zQj) (8) 

VO V 

is the Lienard-Wiechert potential. A gauge transformation has been made so that the field 
is to a good approximation transversal. 

For the elastic form factor F{Q'^) we choose F{Q'^) = exp(Q^i?^), with ^ 



as 
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where 



< > /6 ~ 2.2fm. Alternatively we can set F{Q^) = 1, that is = 0, as the electric 
field outside a spherically symmetric charge distribution is the same as that of a correspond- 
ing point charge. We find numerically the effect of a finite R^ to be rather small, justifying 
this assumption. We still keep it for completeness in the following equations. 

In order to describe the meson production we need an expression for the electromagnetic 
current J {A —>■ A + V). This current can be found by using the vector dominance model 
(VDM), which relates this current to the elastic scattering amplitude V + A^V + A a.s 

J^iQ) = e'^Cvfeii^±,Y) (9) 

with A± = f _L — Q±. ey is the polarisation of the outgoing vector meson, which by assuming 
s-channel helicity conservation (SCHC) is identical to the one of the incomingphoton, see 
[igj ] for details. Cy describes the vector meson content of the photon, see e.g. ^. 
In the following we choose the elastic vector meson scattering amplitude as 

/ez(A^,r) = fo{Y)exp{-AlRl)S{Au.){vu.) (10) 



with Ry ~ 2.2fm to reproduce the slope of the angular distribution. This form agrees also 
with the one proposed in [2^. It has been mainly chosen for ease of calculations, whereas 
the formalism given can be extended to more realistic forms, e.g., based on a full eikonal 
description. A simple extension is possible by using a sum of Gaussians for /g^. 

The imaginary part of /o can be related to the total cross section for vector meson 
scattering on nuclei through the optical theorem. For the energies of RHIC and LHC the 
real part is small (of the order of 10%). While for F 7^ (and asymmetric collisions) 
there is a sensitivity to the phase of /o, for rapidity F = we only need the absolute 
value of |Cy/o|. This we choose in order to reproduce the cross section given in 0, 121I as 
da/dYp{RmC) = 70mb and da /dYj /^{LRC) = 0.75m6. 

Using these expressions we get for the amplitude: 

av(0,f±,F) = -^^--p-exp(-Q;i?^) 



/ 



d^Qi_{Qi_ev, ^2 , n2 



exp(-«Q^6) exp(-i?^({;^ - Q^f) (11) 

with Ql = Ql — Ql = (^)^- Let us first make a short qualitative discussion: if we would 
neglect Q± in exp{—Ry{v± — Q±y), the dependence of \av\^ and d^a/d'^v±dY on v± would 
be of the form exp(— 2i?yf^), which is due to J alone. This coincides with the result for an 
incident photon of zero transverse momentum. The effect of the finite Qj_ distribution of 
the photon is to broaden this distribution. As the width of the Q± distribution depends on 
b via exp{—iQ±b), the effect of this broadening will depend on b. This effect is largest for 
small b, as the perpendicular momentum distribution of the photon is of the order 1/6. 

An analytic approximation for ay can be found in the region of small b, that is if 6 < 1/Qi 
and Qi can be neglected in the photon propagator in Eq. (fTTj). which corresponds to the 
sudden limit. One gets (see [l^ for details): 

ZeCvfoi {b + 2w^Rl)ev 



avib,v±,Y) 



{2nr 

exp{-QfR^) exp{-v\RY 
-(b + 2iv^Rl''^' 

exp ■ 



m'v+R',) ' ^ ■ ^^^^ 

The same final state can be obtained by exchanging the roles of both ions, see Figs, ^a) 
and (b). The corresponding amplitude ay (6, v±, Y) where the photon is emitted from ion 2 



and the "Pomeron" from ion 1 is given by 

a^(6, v^, Y) = j ^^L*(0, Q)Ji{Q), (13) 

where the impact parameter for Af,xt is now 6 = and Uj^ is replaced by m_. The elec- 
tromagnetic current is now for vector meson production on an ion at position h. One 
finds 

= Jf,{Q)exp{-wJ+tAj). (14) 
We find that the exchange amplitude is of the form: 

o^vi^^^-i-^^) = a'vi—b,v±, —Y) exp{—iv±b). (15) 

This has a simple interpretation: In order to exchange the role of the two ions, Y is replaced 
by —Y, the direction of b needs to be reversed and in addition the origin needs to be shifted 
lead., .o .Ke e.t.a pK.e exp,-,...). TM. .e,a«o„ a.o used . QQ. W.. 
ay from Eq. (|TT| we finally get 

a'^\b, v^, Y) = av{b, v±, Y) + e"*^^^ V(-&, -Y). (16) 

The analytic expression in Eq. allows us to discuss some properties of the transverse 
momentum distribution of the process: In the limit v±Ry <^ b one has ay ~ bey and 
tty = —ay, i.e., the amplitudes have a relative sign of —1, leading to destructive interference 
at small b. In the other limit v^Ry ^ b one has ay = ay, i.e., the same relative sign, but 
ay and ay are smaller than in the first case due to the last exponential in Eg. 1121] . The 
transverse momentum distribution is therefore more complex than treated in |irl l22||. 

We can also derive the results starting from the eikonal or Glauber approach to multi- 
photon processes in UPC collisions, see 1^. In this case the scattering amplitude is given 
by 

— * ZTT / — * — * — * 

//i,Glauber(^) = "fc / d^bexp{iKb) {f\ exp{ix{b)) \i) , (17) 

where K = k' — k is the total momentum transfer to the "target" nucleus. The eikonal x{b) 
takes care of all the different elastic and inelastic processes. In our case we have 

X{b) = Xnuc{h) + Xc{h) + Xi{h) + X2{h) + Xv{h). (18) 



The term Xnucib) describes the effect coming from the nuclear interaction between the two 
ions. It can be approximated by exp{ixnuc{bj) ~ Q{b — 2R). The term xc ~ exp(2i?7log(fc6)) 
describes the elastic Coulomb scattering. The last three terms describe the additional elec- 
tromagnetic interactions: the possible excitation of the first and second nucleus and the 
vector meson production. The eikonal phase for the vector meson production process Xv{b) 
can usually be treated in lowest order by expanding the exponential. The second order term 



would describe double p production, which is still sizeable, see 
initial and final states we get 



Bracketing with the 



< f\Gxp{ix{b))\i >^ iexp{ixnuc{b)) exp{+ixcib)) 

< /i|exp(zxi(&))Ki >< /2I exp K2 > 

<V,t2\xvm2> (19) 

with \i >= \ii,i2 > the initial (ground) states of the two ions and \ f >= f2>\V> the 
final states of the ions and the meson. In order for this process to factorise, we made the 
reasonable assumption, that the vector meson production on the excited nucleus is the same 
as the one on the nucleus in the ground state: 

< V,Z2\xv{b)\z2 >^< V,f2\xvmf2 > ■ (20) 

There is a correspondence of these terms to the different semiclassical amplitudes afi{b) 
in Eq. (jSj), which was also explored in [23]. The x is given by 



^^^^^ = / ^^eik^^^' ^' Q^W) (21) 

with ^g-j^ as given in Eq. (jH)) with 6{Qu+) replaced with 5(7(Qo ~ Qz)), which corresponds 
to the expression of the semiclassical amplitude ay in Eq. ((Zj) in the sudden limit. 

The major difference between the two approaches is the presence of the Coulomb eikonal 
Xc(yb). For ^ 1, xc'(&) is a rapidly varying function and one can evaluate the Glauber 
expression by means of the well known saddle point approximation. One obtains the relation 
between the classical impact parameter and the momentum transfer: b = 21]/ K. As the 
momentum transfer to the ion is not measured in our case, one calculates "inclusive" cross 
sections by integrating over K. This gives the same result as in the semiclassical case of 
Eq. (jni) in the sudden limit {Qi = 0), as the Coulomb phase is purely imaginary and the 
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1 as, the approximate analytical 
lo| . They are shown in Fig. |21 



integration over K leads, (see Eq. (10)-(12) of 23]) to a 5-function for the two impact 

parameters in |/^,,Glauberl'- 

We use both the exact expression Eqs. (jH)) and (fTT|). as we 
result Eq. (fT^ to calculate results for the case studied in 
The analytic result is too large in the untagged case, but its shape agrees quite nicely with 
the full calculation. The effect of tagging for small 6 is a shift of the maximum of the curve 
to larger values of v± and a more pronounced interference structure. In Fig. El we also show 
the similar results for J /'^ production at the LHC 
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FIG. 2: The differential cross section da/cPv±dY is shown for production at y = at RHIC 
(Au-Au collisions with 7 = 108). The solid line is the result including the interference, the dashed 
line the result from an incoherent adding of the two processes. The results for the three different 
tagging cases are given in (a), (b) and (c). The approximate result is shown as dotted line only in 
the first (untagged) case. In the other two case it cannot be distinguished from the full calculation. 



Let us summarize our findings: We have put the transverse momentum distribution on 
a firm theoretical basis starting our derivation from the semiclassical approximation or al- 
ternatively from Glauber theory. The meson transverse momentum distribution was derived 
as a function of b and an analytic expression was given. The interference phenomenon was 
derived within this model. As the main outcome we find in this letter that for a good un- 
derstanding of the interference phenomenon a careful study of the transverse momentum 
distribution is essential. Whereas formally the results look similar to the one given in loi|. 
differences appear both in the transverse momentum distribution as a function of b and in 
the form of the interference. This leads to a more complex result in the intermediate v± re- 
gion. This will also be true in the case, where one is moving away from Y = 0. Our findings 
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v_^[MeV/c] v^[MeV/c] v^[MeV/c] 

FIG. 3: The differential cross section da /d?v±dY is shown for J/\& production at y = at LHC 
(Pb-Pb ions with 7 = 3000). The hnes are the same as in Fig. [3 

are important in analyzing the experimental data of STAR and also PHENIX. Results have 
also been given for future LHC measurements, which would be even more interesting for 
J/\E' or even T production. 
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Figure 3: 
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